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Abstract 

Within the spirit of the Hadron Resonance Gas model, we study a representation of the heavy QQ free energy at 
temperatures below the phase transition in terms of the string and heavy-light hadrons. We discuss the string breaking 
phenomenon and the relevance of avoided crossings between the fundamental string and the hadron spectrum. Good 
agreement with lattice data is achieved. 
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1. Introduction 

A fruitful approach to study the confined phase of 
QCD is the Hadron Resonance Gas (HRG) model, in 
which the equation of state is obtained by assuming a 
gas of non-interacting massive stable and point-like par¬ 
ticles, which are taken as the conventional hadrons Ok 
This model has arbitrated the discrepancies between the 
different collaborations in the lattice community 0. 

The effort to study the thermodynamics of QCD has 
led also to important advances in the study of the in¬ 
terquark forces at finite temperature. The vacuum ex¬ 
pectation value of the Polyakov loop is related to the 
propagator of a static quark, and it has been used as 
an order parameter for the confinement/deconfinement 
transition of color charges 0. The free energy of a 
heavy QQ pair in a thermal medium, which can be 
computed from the correlation function of Polyakov 
loops, has important phenomenological consequences 
in the study of quarkonia physics in heavy ion experi¬ 
ments 0]. In line with the HRG approach we have pro¬ 
posed recently Hi a hadronic representation for the 
vacuum expectation value of the Polyakov loop which 
has been quite successful in describing the available lat¬ 
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tice data up to T - 180MeV [5,2]. In this communica¬ 
tion we represent the heavy QQ free energy in terms of 
the HRG, and compare with recent lattice simulations. 


2. Quark potential and string breaking 

One of the most studied quantities in QCD is the 
static energy between heavy sources, such as heavy 
quark systems 181 and also for more exotic states like 
gluelumps | l ). 10], We will present in this section the 
most relevant properties of this interaction which will 
be used later in the study of the heavy QQ free energy. 


2.1. Static energies in QCD 

Inspired by perturbation theory to second order d 
the static interaction between heavy sources A and B in 
QCD is often modeled implementing Casimir scaling. 


- Aa • Ab 



( 1 ) 


where A is the SU(A e ) color group generator corre¬ 
sponding to the representation of the source. The 
Coulomb-like term accounts for the perturbative behav¬ 
ior from one-gluon exchange at short distances, while 
the confining linearly growing term accounts for the 
energy of the string connecting the two heavy quarks. 
Casimir scaling implies a relation between the funda¬ 
mental QQ = 3x3 and adjoint GG = 8x8 color sources 
given by Vngir) = | Vgg(Q- and holds on the lattice for 
N c = 3 [12] and for heavy sources in the fundamental 
representation QQ = 3 x 3 at N c = 3,5,7 lfl3ll . The N c 
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independence of the V q q potential has also been estab¬ 
lished. For N c = 3 one gets 
4a s 

V Q Q(r) = or-—, (2) 

where as — 7r/16 is the effective coupling and <x = 
4k/3 - (0.42 GeV) 2 is the string tension. 


2.2. String breaking 

Using this potential, one can study the transition be¬ 
tween a QQ system and a meson-antimeson system. 
When the separation between the two heavy quarks in¬ 
creases, the energy of the system rises linearly due to 
the string tension term. Eventually this system becomes 
unstable at a critical distance r c when a light qq is cre¬ 
ated from the vacuum, so that a heavy-light meson- 
antimeson MM = ( Qq)(Qq ) channel opens. 

Neglecting the meson-antimeson interaction, the en¬ 
ergies of these channels are 


E Q Q(r ) = m Q + niQ + V Q g(r) , ( 3 ) 

V (Qq)(qQ)^> ~ A e? + \q = 2A » ( 4 ) 

where A Qq = lim^^co - niQj is the mass of the 
heavy-light meson (the mass of the heavy quark it¬ 
self niQ being subtracted). Using these formulas, one 
can estimate the string breaking to be r c - 4Mq/ct ~ 
1.2 findQ, where Mo is the constituent quark mass, 
which is in fair agreement with lattice results [16]. 

More generally, the QQ state can decay into any of 
the many excited states of the meson spectrum after qq 
pair creation from the vacuum, i.e. V;"’"', - (r) = A[V- + 

A^ 1 , or to the baryon spectrum, as long as they have the 
same quantum numbers as the QQ system]^ 



r [fm] 


Figure 1: Adiabatic avoided crossing structure between the QQ fun¬ 
damental state and the MM spectrum. We display as solid lines (blue) 
those states with avoided crossing, while dashed lines (red) stand for 
states with no avoided crossing, i.e. those with different quantum 
numbers than QQ. We have applied a mixing strength of 50MeV. 


3. Heavy QQ free energy 


The free energy is (minus) the maximum work the 
system can perform at fixed temperature. For a QQ pair 
at separation r the free energy is related to the correlator 
between Polyakov loops in the following way [18] 


e -F Q6 (r.T)iT = (TrfQtrlTrfGlO) 1 '). 


(6) 


This is the color averaged free energy, which follows 
from the thermal_average of singlet and adjoint free en¬ 
ergies, see e.g. 1191 1 Several models have been pro¬ 
posed in the literature for the QQ free energy either in 
the confined [2(j r 21 ] or the deconfined phase 22 231. 
We focus on the former case here. 


2.3. Avoided crossings 

In the picture presented above we have assumed a dia- 
batic crossing structure in which the states QQ and MM 
are degenerate at the point r — r c . However, the exis¬ 
tence of a transition potential V between the 
two states, implies the matrix interaction 


V(r) = 


(5) 


v qq(f) Vqq-*mm(f) 

22->mm(G V) 

After diagonalization, the finite energy of the non¬ 
diagonal QQ —> MM interaction lifts the degeneracy, a 
feature called level repulsion, and there emerges the pic¬ 
ture of adiabatic avoided crossing shown in Fig. [T| 11141 


1511 . a phenomenon familiar from molecular physics in 
the Born-Oppenheimer approximation 0- 


3.1. Hadronic representation 

The correlation function in Eq. ([6]» can be written 
as a ratio of partition functions, with and without QQ 
sources placed at a distance r [24]]. As argued in [14], 
in the confined phase such ratio admits an approximated 
spectral decomposition 


e -F QQ (r,T)/T = e -V ae (r)/T + ^ g - 


V^\r)/T 


(7) 


By neglecting the avoided crossing, one finds 


e -F QQ (r,T)/T = e -V QQ (r)!T 




-A JT 


( 8 ) 


2 The mechanism of decaying into baryons implies the creation of _ 

two pairs of light quarks qq, leading to the formation of two heavy- 3 This normalization of the free energy differs by a factor A 2 from 

light baryons with one heavy quark. other conventions found in the literature. 
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where A„ = A ( "i = 


v<«) 


which constitutes a hadronic 


'qQ qQ ’ 

representation of the heavy QQ free energy [14] 15]. In 
Fig. 0A) we display the lattice results of 12511 for the 
heavy QQ free energy along with the calculation from 
Eq. ([8]). We use the spectrum of heavy-light mesons and 
baryons with a charm quark and no strangeness, Nf = 2, 
from the Isgur model lEhll up to A = 3.19 GeV. We 
have taken cr = (0.40 GeV) 2 and as = n/16. We find an 
excellent agreement for the lowest temperatures and for 
all distances. 

The Polyakov loop expectation value can be com¬ 
puted from the large separation limit of the correlator 
between Polyakov loops. From Eq. © one has 


L(T) := lim e ~ F ^ r)KTn = I Y e~^ T . 

r— »oo 7 i 


(9) 


This is precisely the hadronic representation of the 
Polyakov loop introduced in Ref. 151], which was shown 
to provide a good agreement with lattice data for T < 
180MeV. This means that the spectrum from the Isgur 
model already saturates the sum rules at these tempera¬ 
tures. It was shown in [27] that Eq. (|9} can also be ob¬ 
tained in chiral quark models coupled to the Polyakov, 
when one advocates the local and quantum nature of the 
Polyakov loop [28, 29ll . 

A different exercise is to extract from the finite tem¬ 
perature lattice data the zero temperature heavy QQ po¬ 
tential. By using Eqs. ([Sj and (|9}, one gets 


Vqq(Q = -T log 


e -F Q Q(r,T)!T _ e -F eS ^.T)/T 


( 10 ) 


The terms in the argument of the log-function are known 
from lattice data and they are in general T -dependent, 
so that this is a nontrivial check. The result for V 0 g(r) 
is plotted in Fig. [2f B). This quantity turns out to be 
T -independent with a high accuracy, and it is remark¬ 
ably close to the QQ Cornell potential, Eq. Q, giving a 
strong indication of the validity of Eq. ©. 

3.2. Avoided crossing and heavy QQfree energy 


The lattice data of Ref. [251 are computed for 2 light 
flavors with a pion mass of m n - 770 MeV, and as¬ 
sumes the zero temperature QQ potential behavior of 
the free energy at short distances. The recent study of 
Ref. S uses Nf = 2 + 1 with physical quark masses, 
and releases the short distance assumption in the renor¬ 
malization procedure. The lattice results of |[3Cl are dis¬ 
played in Fig. U A) along with our Ay = 2+1 calculation 
assuming the hadronic representation of Eq. ©. The 
most striking observation is that a quite high value of 
the string tension has to be assumed, cr = (0.55 GeV) 2 , 


in order to reproduce the lattice data of 0, while 
cr = (0.40 GeV) 2 worked for the data of [25]. The same 
conclusion is reached by extracting the zero temperature 
QQ potential by using Eq. (ITOb . as shown in Fig.0B). 

In view of this, one may wonder whether this effect 
in the string tension could be effectively due to the ex¬ 
istence of an appreciable mixing between the string and 
the excited states in the lattice calculation of 13011 . To 
investigate this possibility, we take a simple model con¬ 
sisting of only two coupled states, the QQ pair (string 
state) and the lightest meson-antimeson state: 


V(r) = 


4 as 
3 r 

W{r) 


+ err W(r) 
2A 


( 11 ) 


As a first estimate, we assume for the mixing a function 
of the form 

W(r) = Wo e~ mr . (12) 

From a fit of the lowest temperature lattice data of 01 
for the heavy QQ free energy, T = 150 MeV, one gets 


Wo = 1.03(50) GeV, m = 0.85(38) GeV, 
cr = (0.420(10) GeV) 2 , 


(13) 


with^f 2 /dof = 0.036. We have taken 2A = 0.944 GeV 
corresponding to the lowest heavy-light meson state of 
the Isgur model. The result is displayed in Fig. |4]where 
one can see that the mixing mimics a larger value of the 
string tension at distances r < 0.5 fm. The behavior of 
the first excited state at short distances is in qualitative 
agreement with lattice simulations, see e.g. Ref. [16]. 
This result leads to the problem of the identification of 
the parameter m with the mass of some particle. At first 
sight it would be tempting to assume that m should be 
the pion mass, however, the value we get for m is heavier 
by a factor - 2n. 

The model dependence of this result can be relaxed 
by not assuming any functional form for the mixing 
W(r). To do so, for each point of the lattice data |30|] 
at a given temperature we obtain the value of W that 
allows the two states model, Eq. (HU. to reproduce the 
free energy at that separation. The result for W(r) com¬ 
puted in this way is displayed in Fig. [5] A result quite 
similar to Eqs. (TT2l >- ([T3l > is obtained for T < 200 MeV. 
The T -independence of W(r) goes in favor of the relia¬ 
bility of the model of Eq. (fTTV For temperatures above 
150 MeV, one would need to include a huge amount of 
states to reproduce the Polyakov loop. Even after in¬ 
cluding all the available states of the Isgur model, the 
sum rule for the Polyakov loop cannot be saturated for 
temperatures above 180 MeV 11 El. For this reason, 
and to avoid a source of error in W(r) coming from the 
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Figure 2: (A) Color averaged heavy QQ free energy as a function of the separation. The dots are the lattice data for Nf = 2 taken from Ref. f25h . 
The solid lines represent the result by using Eq. {8j with the spectrum of heavy-light mesons and baryons with a charm quark, obtained with the 
Isgur model of Ref. l26ll . (B) Heavy QQ potential at zero temperature as a function of the separation by using Eq. ED and taking the lattice data 
shown in (A) for T / T c = 0.76, 0.81, 0.90 and 0.96. The dashed (black) line is the result of the Cornell potential, Eq. (2}, with cr = (0.40 GeV) 2 . 
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Figure 3: Same as Fig.[2]for the lattice data of Ref. Ell . We 


use Nf = 2+1 for the meson and baryon spectrum and cr = (0.55 GeV) 2 . 


value of F(oo , T), we have used for F(oo, T) the value 
extracted from the lattice data. This is not necessary for 
the lowest temperature T — 150 MeV, as the sum rule is 
saturated in this case. 


4. Discussion and outlook 


We have studied a natural extension of the HRG 
model to the heavy QQ free energy. This quantity has 
been computed in the confined phase of QCD from a 
string and the HRG model with heavy-light hadrons 
and compared with existing lattice computations [25 


13011 giving an excellent agreement for T < 200 MeV. 
The relevance of string breaking at large distances has 
been discussed as well. Finally, we have investigated 


the possible role played by the avoided crossings be¬ 
tween the string and the spectrum of hadrons. 


Our analysis of the avoided crossings leads to pos¬ 
sible contradicting results when analyzing the lattice 
data for the heavy QQ free energy from two different 
groups. While the interaction between the string state 
and meson-antimeson states seems to be negligible in 
the data from 1 251 , a noticeable effect appears in the 
analysis of the more recent data of iBlill . This could be 
well due to the different approaches for quark masses 
assumed in these references, or the different renormal¬ 
ization prescriptions. So, more accurate lattice data and 
better agreement between different groups would be de¬ 
sirable on the one hand. On the other hand, our analysis 
can be improved either by including some subleading 
effects like string excitations or extending the two states 
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Figure 4: With solid lines, string state and first meson-antimeson 
state as a function distance using the model of Eq. m with W(r) = 
Wo e~ mr . The parameters are obtained from a a best fit to the lattice 
data of the heavy QQ free energy of Ref. a at the lowest available 
temperature, T = 150MeV. The dashed lines represent the result 
with this model taking cr = (0.42 GeV) 2 and no mixing, W(r) = 0. 
The dot-dashed line (black) is the result from the Cornell potential, 
V(r) = —j| T + err+ C with cr = (0.55 GeV) 2 and C = -0.3 GeV. 

mixing model with the inclusion of more states in the 
spectrum. Work along these lines is in progress. 
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